IEEE Transactions on Biomedical Engineering

The features of the PPG (photoplethysmography) morphology are known to reflect age-related cardiac and vascular changes. In most contemporary wearables, PPG signals are acquired from distal sites such as the wrist and finger. The superficial temporal artery (STA), accessible at the temple region, is reached via a shorter arterial path from the aortic root than the radial circulation, and may therefore carry hemodynamic and aging information with less distance-dependent attenuation. We hypothesized that the morphology of the PPG at temple region (STA) would show stronger and more numerous age correlates than the PPG at the wrist. To test this, we extracted a common set of 89 pulse-morphology features, spanning raw-waveform timing/amplitude/area measures, ratios among them, derivative-based ratios, and spectral harmonic-ratio features. We compared an in-house temple-worn device which has PPG as one of the sensors, with a publicly available Microsoft Aurora-BP wrist-worn PPG dataset, and tested each feature's association with age. We identified 14 robust age correlates at the temple region, compared to 3 at the wrist. The temple's correlates spanned multiple morphological categories and showed a larger age-association than at the wrist. These results support the hypothesis that the temple region may be a more robust PPG measurement site than the wrist to extract age-related cardiovascular information, which motivates further investigation of temple-based cardiovascular sensing.

Muscle architecture shapes muscle function and changes with age, growth, training and disease, yet quantifying three-dimensional (3D) muscle architecture in vivo remains challenging. We introduce a hybrid fascicle tractography approach for freehand 3D ultrasound data that accurately reconstructs 3D muscle fascicles with respect to an objective, anatomically relevant coordinate system defined by the muscle's central aponeurosis. The hybrid approach combines Hessian-based fascicle detection with wavelet-based refinement to generate volumetric fascicle orientations. In a synthetic dataset with known ground truth, fascicle orientations and lengths were estimated with errors of [≤]2{degrees} and ~1.5%, respectively. In vivo, the approach detected physiologically plausible fascicle lengthening in the human tibialis anterior following a passive plantar flexion rotation, whereas diffusion tensor imaging of the same muscle did not. The proposed method enables anatomically relevant, objective and non-invasive quantification of 3D muscle architecture in vivo, providing a practical framework for applications in clinical and applied muscle physiology.

Surface electromyography (sEMG) offers a silent and wearable input modality, but its practical use is limited by variability across users and recording sessions. This study presents a compact CNN- Transformer model for decoding isolated handwritten digits from eight-channel sEMG signals. The model combines trainable signal preprocessing, convolutional feature extraction, and Transformerbased temporal modeling. It was evaluated on ten recordings from five participants using recordingseen classification, leave-one-recording-out (LORO) generalization, and few-shot adaptation. The model achieved a mean macro F1 score of 0.924 {+/-} 0.059 in the recording-seen setting and 0.619 {+/-} 0.252 under zero-shot LORO evaluation. Adaptation using two labeled trials per digit increased macro F1 to 0.828 {+/-} 0.112, while ten trials per digit achieved 0.925 {+/-} 0.053. The proposed architecture also outperformed classical and neural baselines in the controlled LORO benchmark. These results indicate that compact CNN-Transformer models, combined with lightweight target-recording calibration, provide a promising basis for adaptive sEMG-based input systems.

Neural temporal point processes (NTPPs) are powerful tools for modeling sequences of timestamped events with statistical temporal structure. Density-based NTPPs, in particular, are an interesting opportunity to merge the universal function approximation capability of neural networks with a defined statistical model in a way that has many potential applications. We demonstrate one such application to heartbeat dynamics, a physiologic point process. We specifically apply a lognormal mixture NTPP to compute instantaneous estimates of the mean and standard deviation of beat-to-beat intervals. We compare our results to the state of art (Barbieri et al.) point process model for heartbeat dynamics, which uses a more physiologically rigorous inverse Gaussian model. We find that the NTPP model maintains reasonable accuracy while improving upon robustness to noise.

Cerebral hemodynamics are difficult to monitor continuously outside the laboratory. Optical head-worn wearables have been proposed for tracking cerebral blood-flow signals, but they require comparison with an established cerebrovascular reference before they can be interpreted. We evaluated a temple-worn optical wearable, Temple, that outputs a proprietary, dimensionless Brain Flow index, intended as a proxy for relative changes in cerebral hemodynamics, against transcranial Doppler (TCD) ultrasound, which measures blood-flow velocity in the middle cerebral artery (MCAv). Twenty-three healthy adults completed two physiological challenges that elicit distinct and acute cerebral hemodynamic responses: a cycle-ergometer exercise protocol and a stand-to-supine postural transition protocol. Twenty participants were analyzed per protocol. The Brain Flow index tracked MCAv in both protocols, with significant within-subject temporal correlations (median Pearson r = 0.795 and 0.799 for exercise and postural transition; p < 0.001) and directionally concordant, statistically significant transition responses for both increases and decreases in flow. Bland-Altman analysis of the normalized transition responses showed small mean biases between the two devices, consistent with similar relative response shapes. Because both signals were standardized within session before this comparison, it addresses the shape of the relative change rather than agreement in absolute units. The Brain Flow index reproduced the direction and time course of MCAv under both perturbations, including the postural transition, where heart rate moved in the opposite direction. Further studies using complementary modalities and additional cerebrovascular reactivity challenges are required to establish clinical use cases and cerebral specificity of the Brain Flow index.

Pathological tremor is a neurological condition that impairs fine motor tasks, affecting 1% of the general population and 4% of the elderly. Tremors arise when muscles micro-oscillations synchronize and phase lock, typically within a 4-12 Hz frequency range. Administering beta-blockers can reduce tremor severity, but doses are hard to personalize, with heavy doses of propranolol correlating with low blood pressure, dizziness, and nausea. In this project, we aimed to model tremor and create a closed-loop control framework to suppress tremor amplitude while minimizing pharmacological dependence. Because side effects constrain the use of pharmacological suppression alone, we investigated noninvasive neuromodulation. We used vibrotactile stimulation (VTS) to disrupt pathological tremor synchronization and reduce oscillatory amplitude. We hypothesized that tremor suppression involving VTS followed a nonmonotonic relationship, tested by determining whether maximum relief requires an adaptable framework. The procedure consisted of constructing a propranolol-reduction simulation by implementing a Hill curve, where we calculated and utilized tremor reduction, heart rate (HR) drop, and blood pressure (BP) drop. We then built a device to capture tremor-related data and create vibration using two linear resonant actuator (LRA) coin motors. We connected it to a microcontroller, where we determined optimal vibration frequencies through a feedback loop. Across 50 trials, VTS alone reduced tremor amplitude by an average of 37.3%, reducing the propranolol dose needed to reach 50% total tremor reduction by 71.9%, lowering the modeled blood pressure drop from 38.1 to 18.9 mmHg. This device demonstrates proof-of-concept for a nonmonotonic tremor-vibration relationship to reduce dependency on propranolol in the treatment of pathological tremor. These propranolol dose-reduction estimates are derived from computational simulation and have not been clinically validated; they are not intended as a recommendation to alter prescribed medication.

Breast cancer survivors frequently experience upper-limb impairments, making continuous monitoring essential for effective rehabilitation. We propose REINA (Recognize-Then-Infer Wearable-to-App AI Framework), a two-stage deep-learning approach for remote monitoring of motor function during breast cancer rehabilitation using wearable-device data. Inertial measurement unit (IMU) signals from wearable devices are first used to recognize physical activities via supervised learning, followed by an activity-specific recurrent neural network (RNN) to infer corresponding electromyography (EMG) signals. REINA establishes reliable inference of neuromuscular activity from wearable IMU data, enabling real-time, cost-effective assessment of motor function recovery in real-world settings.

Objective: Low-intensity transcranial magnetic stimulation (LI-TMS) is being investigated as a gel-free alternative to transcranial electrical stimulation (tES), but existing systems remain almost exclusively single-channel and cannot electronically steer the induced electric field. We present the design, modeling, and experimental measurement of a wearable whole-head, multichannel, steerable LI-TMS array. Methods: The system comprises a 102-channel conformal coil array with independently controlled drivers capable of arbitrary waveform synthesis, together with a boundary element fast multipole method (BEM-FMM) framework that computes the coil currents required to produce prescribed cortical field patterns. A 12-channel prototype was characterized by coil-current, electric-field, and thermal measurements. Results: The prototype produced a peak primary electric field of approximately 1 V/m measured in air 4 cm from the inner helmet surface. Whole-array modeling attained cortical fields of up to 1.5 V/m, reproduced the field distribution of a clinically validated low-intensity stimulator to within 3%-5%, and demonstrated focal targeting of the dorsolateral prefrontal cortex, simultaneous delivery of electric field to the default mode network nodes, and synthesis of electric fields following the traveling alpha wave. Conclusion: Electronically steerable, whole-head LI-TMS is feasible using accessible microprocessor-controlled power electronics. Significance: The array reaches the cortical field regime of tES without scalp contact or the associated shunting of current through the scalp, offering a route to testing network-level, phaselocked weak-field neuromodulation.

Wearable magnetoencephalography (MEG) enabled by optically pumped magnetometers (OPMs) promises improved comfort and motion tolerance. This is particularly beneficial when measuring brain activity in children who cannot sit still for long periods of time. Compared to cryogenic MEG, wearable MEG allows larger head movements, but they result in artifacts due to uncompensated background fields and reduce source localization accuracy. Spatial filtering methods can partially compensate these motion-induced artifacts, but they are most effective when used in combination with background field nulling. This is because accurate spatial filtering relies on an accurate estimate of the sensor gain and orientation of its sensitive axis. Through simulations, we first deduce the target residual background field that is necessary for accurate dipole localization (< 1 cm) in the presence of head movements. Using our open-source printed circuit board (PCB) coils, we develop a method to dynamically null the background field. We demonstrate that our dynamic field nulling method allows improved localization of somatosensory evoked fields (SEFs) by maintaining the background field below the target residual fields established in the simulations. Our study highlights the importance of tracking both the background field and the head position relative to the background field for quality assurance in wearable MEG.

Coronary flow waveforms have a distinct diastolic-dominant shape with periods of low or retrograde flow during systole. While the general waveform shape has been attributed to complex interactions between cardiac and vascular mechanics, there is limited research into the variability in coronary flow waveforms and what this variability may reveal about cardiac function. This work presents a shape analysis of left anterior descending artery (LAD) flow waveforms using Fourier transforms and Singular Value Decomposition (SVD) performed on baseline data collected from 32 pigs. Pigs included in the study reflect two breeds (Ossabaw and Yorkshire) and three different experimental conditions (lean-control, lean-paced, and obese-paced). Fourier transforms were used to decompose the waveforms into 15 harmonics for each pig. An SVD analysis is then used to extract temporal patterns of the waveforms. Correlations between pig-specific coefficients for the SVD modes and clinical metrics were used to investigate physiological explanations of LAD waveform variability. Temporal LAD flow patterns of the second SVD mode are significantly correlated with heart rate. The third SVD mode significantly correlates with mean blood pressure and maximum hyperemic flow. Furthermore, the fourth SVD mode is weakly correlated with left-ventricular end diastolic pressure and endocardial-epicardial flow ratios. This work demonstrates that LAD flow waveforms can be broken down into temporal patterns that correlate with physiological features. Furthermore, this shape-analysis method allows for waveform reconstruction and simplifies visualization of the temporal patterns identified using SVD, an advantage over existing methods that focus on characterizing flow waveforms by points of interest.

Physical human-device interaction is central to many emerging technologies in neurorehabilitation and assistive robotics, but simulation-based research in this area remains fragmented across musculoskeletal models, assistive-device representations, task definitions, and controller-development workflows. This fragmentation limits the accessibility, reproducibility, and extensibility of studies on prostheses, exoskeletons, wearable rehabilitation devices, and related human-device systems. Here we introduce MyoAssist 1.0, an open-source framework for neuromechanical simulation of physical human-device interaction built within the MyoSuite ecosystem. MyoAssist organizes each simulation environment as a composed human-device-task system that combines compatible musculoskeletal, assistive-device, and task-scenario components through a shared composition pipeline. The current release includes 15 assistive-device models spanning gait assistance, upper-body support, manipulation, and seated mobility and supports compatible musculoskeletal models ranging from reduced lower-limb models to a 416-muscle full-body model. These human-device systems can be simulated within the broad task scenarios provided by MyoSuite, while MyoAssist adds locomotion-specific task scenarios with configurable terrain and target-velocity conditions for gait-assistive studies. MyoAssist also provides two complementary controller-development frameworks: a reinforcement-learning framework for training adaptive policies and a controller-optimization framework for tuning structured, interpretable human and device controllers. Both frameworks operate on the same simulation environments and provide standardized evaluation outputs for inspecting, comparing, reusing, and extending learned and structured control strategies. By integrating modular human models, assistive-device models, task scenarios, and training workflows under a shared open-source interface, MyoAssist aims to lower the barrier to reproducible simulation-based research and to support collaborative development of assistive technologies for neurorehabilitation and physical human-device interaction.

Purpose: Patient identity management is fundamental to healthcare information systems, as identification inconsistencies can compromise patient safety, data integrity, and clinical workflow efficiency. Reliable linkage of medical images acquired across different imaging modalities remains challenging because of variations in image appearance, acquisition geometry, and imaging characteristics. In this study, we developed an automated patient identity verification framework for multimodal medical imaging using deep metric learning and Data-Augmented Domain Adaptation (DADA). Methods: The proposed framework learned modality-invariant patient representations from labeled source-domain data while leveraging unlabeled target-domain data to mitigate cross-modality distribution shifts. Chest radiographs and computed tomography (CT) scout images obtained under routine clinical conditions were retrospectively collected and used for evaluation. Verification performance was assessed using receiver operating characteristic (ROC) analysis, with the area under the ROC curve (AUC) used as the primary performance metric. Results: The proposed framework achieved consistently high verification performance across all evaluation conditions, with AUC values ranging from 0.9997 to 0.9998. Similarity-score distributions demonstrated distinct separation between same-patient and different-patient image pairs despite substantial differences between imaging modalities. Conclusion: These findings indicate that patient-specific anatomical representations can be preserved across heterogeneous imaging domains through metric learning and domain adaptation. The proposed framework may serve as a practical infrastructure component for patient identity management, multimodal data integration, quality assurance, and patient safety applications within healthcare information systems.

Surface electromyography (sEMG) is a robust non-invasive modality for human-machine interaction, yet its application remains largely limited to coarse motor tasks such as grasping or rotation. The decoding of fine motor skills, specifically handwriting, remains a challenging problem with potential relevance for prosthetic control and natural communication interfaces. In this work, we explore a Transformer-based alternative to classical signal-processing pipelines that treats multi-channel sEMG signals as complex time series. We introduce DualMyo, a specialized model integrating Patch Embeddings and Rotary Positional Embeddings (RoPE) to capture the intricate spatio-temporal dynamics of myoelectric activity. Our experimental results show strong intra-session performance. Furthermore, we address the inherent challenges of signal drift and sensor displacement in cross-session applications. We show that a lightweight fine-tuning strategy of 10 epochs enables DualMyo to effectively adapt to session variability, achieving approximately 91\% accuracy with two examples per digit. These findings provide a promising step toward adaptive sEMG-based handwriting interfaces, although further validation is required for real-time and multi-subject deployment and neuromuscular control.

Cerebrovascular pulsatility drives measurable brain tissue deformation and has been associated with ageing and a range of neurological disorders. Yet how pulsatile haemodynamic forces are transmitted through deformable cerebral arteries into the surrounding brain remains poorly understood, particularly in anatomically realistic vascular geometries. Existing computational approaches have largely treated cerebral fluid and tissue mechanics separately or relied on idealised geometries, limiting our ability to determine how vascular anatomy simultaneously governs intraluminal haemodynamics and extravascular mechanical loading. Here, we develop an image-derived three-dimensional computational framework that jointly resolves pulsatile blood flow, arterial wall deformation and surrounding brain tissue motion in representative cerebral arteries. Four arterial segments, including the middle cerebral artery, middle cerebral artery bifurcation, basilar artery and internal carotid artery, are reconstructed from high-field (5 Tesla) magnetic resonance imaging data of a healthy subject. A finite-deformation fluid-structure interaction model is established by coupling non-Newtonian blood flow, hyperelastic arterial wall and hyper-viscoelastic brain tissue. The predicted tissue response is benchmarked against in vivo magnetic resonance elastography measurements of cardiac-induced volumetric strain over a cardiac cycle. Results reveal spatially localised arterial and tissue deformation whose magnitude and distribution are strongly governed by vascular geometry and wall thickness. Among the segments examined, the internal carotid artery exhibits the largest deformation response, while reduced wall thickness increases strain transmission into the surrounding tissue. Geometrically complex regions also exhibit greater spatial heterogeneity in near-wall haemodynamic metrics. These findings demonstrate that cerebral vascular anatomy simultaneously shapes intraluminal haemodynamics and extravascular mechanical loading. By integrating image-derived vascular anatomy, coupled blood-vessel-brain mechanics and in vivo benchmarking within a unified framework, this study provides a mechanically consistent reference for healthy cerebral pulsatility and establishes a foundation for quantifying how blood-vessel-brain interactions are altered under pathological conditions.

Static force plates (FP) are the gold standard for measuring ground reaction forces (GRF) and computing joint moments through inverse dynamics in gait analysis. However, they are restricted to controlled environments, and the number of steps analyzed is limited by the plates embedded in the floor. To address these limitations, portable solutions such as sensorized insoles, socks, or shoes have emerged. Yet, creating wearable systems capable of measuring three-dimensional GRF in real-world conditions remains challenging. Current sensorized shoes often incorporate thick sensors (up to 2 cm), reducing usability and limiting their application in pathological populations or dynamic tasks like running. This study evaluates the usability of ShokacShoes, a novel sensorized shoe integrating three thin, three-dimensional force sensors, and explores its potential as a Wearable Force Plate (WFP). Eight healthy participants performed slow, natural, and fast walking using two insole configurations. Force and temporal metrics were derived from WFP and FP data. Results indicate that WFP enables accurate step segmentation and detects significant effects of speed and insole type on temporal and force metrics, confirming its reliability under different walking conditions. Comparisons with FP revealed differences in force metrics and signal morphology, though temporal parameters remained consistent. These results are likely due to sensor quantity and positioning. Thereby, ShokacShoes represent a valid solution capable of measuring three-dimensional forces within commercial footwear. Future work will focus on validating the applicability of a new version of ShokacShoes against gold-standard FP in a comprehensive validation study involving diverse real-world scenarios and pathological conditions.

Objective: Accurate volumetric analysis of the brain and cerebrospinal fluid (CSF) is essential for monitoring hydrocephalus, a significant pediatric neurological condition. While computed tomography (CT) provides high-quality volumetric assessment, its associated ionizing radiation poses risks, especially for children. Low-field magnetic resonance imaging (LF-MRI) offers a safer and more accessible alternative, particularly in resource-constrained settings. However, its lower resolution and increased susceptibility to structural distortions make accurate segmentation challenging. This study aims to demonstrate that reliable volumetric measurements can be obtained from LF-MRI that are comparable to CT, enabling safer and more frequent monitoring of infants with hydrocephalus. Approach: We propose EnSegNet-Cross, a cross-modality, enhancement-aware segmentation network for brain volume analysis using LF-MRI. The framework leverages high-fidelity CT data during training but requires only LF-MRI during inference. At the core of the framework is a novel cross-modal topological penalty designed to minimize discrepancies between predicted LF-MRI and CT structures. A central contribution is the integration of a three-dimensional topological loss based on persistent homology, which penalizes topological discrepancies in CSF regions, specifically CSF holes formed by enclosed brain parenchyma, between CT and LF-MRI segmentations. Incorporating these structural priors facilitates generalization across heterogeneous clinical cases while eliminating the need for CT data during inference, resulting in more anatomically coherent and topologically faithful segmentations. Main Results: On a curated cohort of infants with hydrocephalus who had paired LF-MRI and CT scans, including cases with infectious and non-infectious causes, EnSegNet-Cross consistently outperformed state-of-the-art machine learning alternatives. It achieved the highest Dice score of 0.8532 plus/minus 0.03 and Volume Score of 0.9318 plus/minus 0.03. The method also demonstrated robust performance in challenging cases with confounding factors, achieving a Dice score of 0.8340 plus/minus 0.03 and a Volume Score of 0.9111 plus/minus 0.05. By leveraging CT-derived topological priors, EnSegNet-Cross successfully handled anatomically complex scenarios in which conventional models failed. Significance: EnSegNet-Cross provides a reliable and interpretable solution for brain and CSF segmentation, particularly in complex cases of hydrocephalus. This study demonstrates that high-fidelity volumetric estimates can be achieved using only LF-MRI, facilitating frequent, radiation-free monitoring. By bridging the fidelity gap between low-quality LF-MRI and high-resolution CT through clinically grounded enhancement and topological supervision, EnSegNet-Cross offers a robust clinical tool for brain volumetric analysis in infants with hydrocephalus using LF-MRI.

Continuous wave diffuse correlation spectroscopy (cw-DCS) is a noninvasive optical technique to monitor the tissues blood flow changes. This technique measures the tissue blood flow index (BFI) by evaluating the decay rate of the autocorrelation function. The derived BFI is proportional to mean squared displacements of the red blood cells considered as the fast-dynamic scatterer component of tissue in time. However, biological tissue contains static scatterer component and slow-dynamic scatterer component which affect the decay rate of autocorrelation function and as a result the derived BFI. In this study, we assessed the fractional contribution of static, slow-dynamic and fast-dynamic scatterer components of a medium in the flow index derived by cw-DCS. The measurements performed on Agar-based phantom with tube showed that presence of static scatterer component and slow-dynamic scatterer component led to substantial underestimation ({approx} 123%) of the flow index derived by Siegert relation, compared to effective diffusion coefficient of fast-dynamic scatterers components derived by modified Siegert relation and bi-exponential model. The less underestimation was observed for the corresponding parameters obtained from the liquid phantom measurements ({approx} 25%) as well as during the forearm occlusion test and respiratory challenges ({approx} 16% - 26%).

People with lower limb amputation lack somatosensory feedback from their prosthesis, and this loss contributes to functional deficits, including balance and gait impairments. Recent advances in neuroprostheses have demonstrated that electrical stimulation of sensory nerves in the residual limb and spinal cord can restore lost sensations in the lower limb. To maximize the efficacy of these somatosensory neuroprostheses, the restored sensations should be intuitive, seamlessly integrating into the sensorimotor network. However, it is challenging to quantify the intuitiveness of these evoked sensations. Recent studies have proposed using crossmodal congruency effect (CCE) tasks for this purpose in people with upper-limb amputation. The current study tests the feasibility of the CCE task for assessing the intuitiveness of sensory feedback in the lower limb. We hypothesized that CCE score would reliably differentiate between a more natural (pneumatic) sensation and a less natural (electric) sensation at two locations: the knee and the foot. Across fifteen able-bodied individuals, we observed that the CCE task differentiates sensory modalities at the knee, but not at the foot. Identification of external factors affecting the CCE is needed before it can be implemented to measure intuitiveness of sensory feedback in lower-limb amputees.

Transcranial magnetic stimulation (TMS) excites neurons noninvasively by electromagnetic induction and is used in neurophysiology research and in approved therapy for depression. Commercial stimulators cost tens of thousands of dollars. Existing open-source designs are either low-energy and unvalidated or rely on expensive switches and laboratory infrastructure. We present a monophasic, fixed-pulse-shape TMS device built at a parts cost of ~USD 700 which, under specific modeling assumptions, can exceed average human motor thresholds. Our design assumes access to basic, off-the-shelf equipment such as a 24 V power supply unit, an oscilloscope, and a few basic tools. The device charges a 230 F film-capacitor bank and discharges it through a self-wound figure-of-eight coil using a thyristor, producing a fixed pulse with a positive lobe lasting approximately 90 s. A Zero-Voltage Switching (ZVS) driver-based charging circuit charges the capacitor bank up to 1460 V from a 24 V bench supply. Three galvanically isolated voltage domains, redundant interlocks, and passive and active discharge paths help mitigate the safety risks involved with handling lethal energy levels. We also present a low-cost way to characterize the device by reconstructing coil di/dt from pickup-coil dB/dt maps to estimate the induced cortical E-fields. At the maximum capacitor voltage, the recovered maximal di/dt is 110.86 A/s, giving estimated 99.9th percentile cortical E-fields of 159 V/m at Oz and 196 V/m at C3 on an example anatomy. Although not yet approved for clinical trials and routine stimulation, the device demonstrated the possibility of a cost-effective TMS unit.

Abstract Optically pumped magnetometer (OPM) sensors can be placed closer to the scalp than conventional superconducting quantum interference devices (SQUID), resulting in larger magnetoencephalography (MEG) signals from neuronal activity. For event-related sensor data, such as epileptogenic activity or sensory and motor evoked responses, however, OPMs and SQUIDs often differ less in signal-to-noise ratio (SNR) than in signal magnitude. We examined two factors contributing to the relative SNR: the dependence of the signal magnitude on source depth and the effect of scalp-to-sensor distance on the noise level. Simulated MEG data for a current dipole in a spherical head model confirmed that on-scalp sensor placement delivers the largest SNR gain for superficial sources. Depending on the relative overall noise level, there may be a crossover source depth at which SNR is equal for on-scalp and off-scalp sensors and beyond which off-scalp sensors achieve higher SNR. Analysis of the equal-SNR source depth in different-sized spherical head models indicated that, for a given relative noise level, the proportion of the brain where SNR is higher in OPM than in SQUID was larger in small head models, supporting the benefits of OPMs in pediatric studies. To experimentally evaluate noise contributions of brain and non-brain origin to the SNR, we recorded somatosensory evoked fields (SEFs) at varying scalp-to-sensor distances. Generally, both the evoked response magnitude and the noise level were lower when the sensors were further away from the scalp; consequently, the SNR depended less than the signal magnitude on the scalp-to-sensor distance. Comparison of power spectral densities (PSDs) at different sensor-to-scalp distances allowed us to identify whether the dominant noise source was of brain or non-brain origin at different frequency bands. Overall, the results highlight complementary properties of OPMs vs. SQUIDs in terms of SNR, which is of interest when optimizing MEG experiments for specific subject populations and brain regions.